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Abstract 

We derive semi-analytic expressions for the analytic continuation of the Mellin transforms 
of the heavy flavor QCD coefficient functions for neutral current deep inelastic scattering 
in leading and next-to-leading order to complex values of the Mellin variable N. These 
representations are used in Mellin-space QCD evolution programs to provide fast evalua- 
tions of the heavy flavor contributions to the structure functions F2(x, Q 2 ), Fl(%, Q 2 ) and 
9i(x,Q 2 ). 



1 On leave of absence from IHEP, RU-142281 Protvino, Russia 



1 Introduction 



A major goal in investigating the scaling violations of the structure functions as measured in 
deeply inelastic scattering is to determine the QCD scale parameter Aq CD , or equivalently the 
strong coupling constant a s (Q 2 ), at a typical reference scale Q 2 . The QCD renormalization group 
equations for mass factorization at space-like momentum transfer are solved referring to either 
rr-space [l] 2 or Mellin-space [3,4] representations. 3 The formulation of the evolution equations in 
Mellin space has the advantage that only ordinary differential equations, unlike integro-differential 
equations in x-space, have to be solved. Moreover, to maintain factorization-scheme invariance 
the evolution equations have to be consistently expanded in the strong coupling constant, which 
allows an analytic integration of the differential equations, i.e. a complete analytic solution. The 
only numerical element consists in the inverse Mellin transform of the solution to rr-space via a 
contour integral around the singularities of the solution in the complex iV-plane. The solution 
in fixed order perturbation theory has isolated poles on the real axis bounded by some positive 
number from above. Resummations of small-x contributions can be easily implemented in this 
approach in a consistent way [6]. In the QCD fit only the parameters of the non-perturbative 
input distributions and Aqcd have to be varied and the evolution kernels consisting out of 
the anomalous dimensions and the Wilson coefficients to the respective order in the coupling 
constant need not to be recalculated but are kept fixed as (tabulated) Q 2 -dependent arrays along 
the contour used for the numerical Mellin inversion. This makes QCD fits using Mellin-space 
programs particularly fast. 

The analytic solution of the evolution equations in Mellin space allows to study factorization- 
scheme invariant evolution equations of different kind [4,7]. These equations describe the evo- 
lution of physical observables. In the non-singlet case a flavor non-singlet combination of struc- 
ture functions forms the input distribution, which can be measured at a starting scale Qq. 
The usual flavor singlet evolution can be mapped into the evolution of two structure func- 
tions as e.g. F 2 (x,Q 2 ) and F L (x,Q 2 ), F 2 (x,Q 2 ) and dF 2 (x, Q 2 )/dln(Q 2 ), or gi(x,Q 2 ) and 
dgi(x, Q 2 )/d\n(Q 2 ), [4]. Here the non-perturbative input distributions are formed by the mea- 
sured structure functions at a scale Qq, and have therefore not to be determined in the QCD- 
analysis. The evolution kernels in x-space are complicated multiple Mellin convolutions of Wilson 
coefficients, splitting functions P^\z) and of some Mellin-inverse splitting functions and Wilson 
coefficients in general, see [8]. These physical evolution kernels can be easily calculated in analytic 
form in Mellin-iV space unlike in x-space, where, as is well known, already the Mellin-inverse 
of leading order splitting functions takes a complicated form [9] and the convolutions to be per- 
formed were never done analytically. Due to the sizeable contribution of charm to the structure 
functions F 2 (x, Q 2 ) and Fl(x, Q 2 ) in the kinematic regime of HERA (x > 10~ 4 , Q 2 > 5 GeV 2 ) an 
exact treatment requires the account for heavy quark effects on F 2 (x, Q 2 ) and Fl(x, Q 2 ) which is 
easiest described in analytic form in Mellin-iV space. The standard QCD analysis and fits using 
scheme-invariant evolution provide different methods to determine the strong coupling constant 
cc s (M§). In the upcoming high-precision analyses, which aim on measurement errors of 0(1%), 
it is of particular importance to diminish all theoretical and conceptual uncertainties. For this 
purpose the use of both these methods will be essential. 

So far no specific representation of the heavy flavor Wilson coefficients in Mellin-iV space, 
in the way used for those of the light flavors, were derived. It is the aim of the present paper 
to present these parameterizations in terms of a semi-analytic approach. The representation is 

2 A series of these codes uses orthogonal polynomials to diagonalize the evolution equations [2]. 
3 For numerical comparisons of different evolution codes see [5]. 
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based on meromorphic functions in the Mellin variable N. The £ = Q 2 /m 2 dependence is given 
by a series of numerical expansion coefficients until a desired accuracy is reached. Interpolations 
in the variable £ are used along the inversion contour to keep the analysis as fast as possible. The 
present parameterizations enable to perform scheme-invariant evolutions including the effect of 
heavy flavor coefficient functions. 

The paper is organized as follows. In section 2 we summarize the basic notations to describe 
the neutral-current heavy flavor contributions to deep-inelastic scattering up to next-to-leading 
order [10-14]. In section 3 we present the parameterization of the different contributions to the 
QCD Wilson coefficients for heavy flavor production for complex values of Mellin-iV and apply 
these representations to calculate the heavy flavor structure functions using sample densities for 
the light partons, including the study of their accuracy. Section 4 contains the conclusions. 



2 The Heavy Flavor Structure Functions 

The Wilson-coefficients for heavy flavor production in deeply inelastic scattering depend on two 
different scales, which can be represented choosing different kinematic variables. Here we follow, 
for convenience, the notation of Ref. [14] and use the variables £ and rj : 

m A 4ra 2 

s denotes the cms energy squared of the heavy quark system, Q 2 = — q 2 the four-momentum 
transfer squared, and m the heavy quark mass. The momentum fraction z carried by the struck 
parton in the nucleon is 



Q 2 . C/4 



Q 2 + s 1 + 7] + f /4 ' 
and the cms velocity of the heavy quarks is 



z e 



Q 

X, 



Q 2 + Am 2 



(2) 



/ 4 771^ Z / Tj 

The heavy flavor structure functions are given in leading (LO) and next-to-leading order (NLO) 
by 

F k (x,Q 2 ,m 2 ) = F k LO (x,Q 2 ,m 2 ) + F^°(x,Q 2 ,m 2 ) , (4) 
G k (x,Q 2 ,m 2 ) = G L k °(x,Q 2 ,m 2 ) , (5) 

with k = 2, L for the unpolarized structure functions and k — 1 for the polarized structure 
functions. In the latter case the NLO corrections are not yet calculated. The heavy flavor 
contributions to the structure function g2(x, Q 2 ) are obtained at leading twist from G\(x, Q 2 , m 2 ) 
via the Wandzura-Wilczek relation [15]. The leading order contributions are 

H™(x,Q 2 ,m 2 ) = -^UoxVq P X - K (-, A , (6) 

TTTfl Jx Z \Z / 

with H — F,G, h — f,g, f g (z,fi 2 ) = G(z, /i 2 ), g g (z, fi 2 ) = AG(z,/i 2 ) the unpolarized and polar- 
ized gluon distribution functions, resp., and an implicit z dependence of the Wilson-coefficients 
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c^jj (r),£). eg denotes the electric charge of the heavy quark, a s = a s /(47r) the strong coupling 
constant, and fj, 2 the factorization scale. The LO Wilson coefficients are [10-12] : 



CM 



7T, 



ix 1 

2 / (TT^+W 



f 2 

\2 , S 



{2[r ] (l + r ] )f 2 -L(r ] )} , 



(7) 



1/2 



(l + 77-e/4) 2 + l+77 



1 + V 



+ 



2(1 +r]Y + ^- + l + 2r] 



4vr 1 f |~ 



3(1 + ^7) -7 



£ / l + r ? + e/4 iVl + 77 
with T/ = 1/2 in SU(N) and 

[v} ~ [(I + v) 1/2 - v 1/2 . 

The NLO contributions F^ LO (x,Q 2 ,m 2 ) are given by : 



(8) 
(9) 
(10) 



'11) 



F^°(x,Q 2 ,m 2 ) = 



nm z Jx 



2max dz 

x z 
3 



e lfg ( 



cS(e,^)+C(e^)ln 

2 \ 



to" 



cS(e,^) + cg(e^)ln 

4^)+4!^) 



to" 



(12) 



with d^ q (^,r]) = 0. 4 In [14] the functions cj^j, d^] and c^j,d^j are parameterized in terms of 
analytic expressions for the threshold and asymptotic region and a tabulated numerical term 
interpolating between both. 



3 Parameterization 

The parameterization of the Mellin transform for the Wilson coefficients in section 2 shall be 
carried out for complex values of the Mellin variable N. Complete analytic results for the 
Mellin transforms of the heavy flavor coefficient functions are difficult to obtain even in lowest 
order in QCD 5 . A very efficient method in representing higher order functions to high preci- 
sion, cf. [16], consists in the MINIMAX-method, see also [17]. This polynomial representation 
uses Chebyshev-polynomials in the approximation. If compared to the Taylor-expansion of a 
function the coefficients are adapted in such a way that the convergence is significantly better 
comparing polynomials of the same degree. This method has been also successfully being used in 
representing analytic continuations [18] of Mellin transforms occurring in massless higher order 
calculations [19] before. 

4 Thcsc functions emerge in case one would like to include photo-production as well. 

5 As shown in Ref. [11] for the unpolarized leading order coefficient functions analytic results for the Mellin 
transform can be calculated for integer N are given in a form which does not allow analytic continuation. 
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The heavy quark structure function (4) are Mellin convolutions 

F(x) — (A<S> C)(x) — ( dxi [ dx 2 A(x 1 )C(x 2 )8(x - x x x 2 ) , (13) 

Jo Jo 

where the parton densities have support X\ e [0, 1] and the coefficient function x 2 e [0,£/(£ + 4)]. 
(13) can be diagonalized by the Mellin transform 

M[F(x)](N) = C dxx N - 1 F(x) (14) 
Jo 

M[F(x)](N) = M[A(x 1 )](N)-M[C(x 2 )](N) , (15) 

accounting for the support in C{x 2 ). The Wilson-coefficients C(z,£) are now parameterized as 
polynomials in z in the range [0,p] with p — + 4) at fixed values of £ using the MINIMAX- 
method to obtain optimal convergence. Sometimes it is useful to study the function 

C MINIMAX (z,0(p- zY = J2a k (p)z k (16) 

k=0 

instead of C(z,£) to improve the convergence. 

The Mellin transform of C(z,£) is then obtained by a Riemann-Liouville fractional inte- 
gral [20] 

M[C{z^)]{N) = j^a k {p)p N+k - K B{N + kA-K) . (17) 

k=0 

Here, B(a, b) denotes the Euler Beta-function. As evident, (17) defines a meromorphic function 
in N with poles on the real axis bounded from above. For k — 0, B(N + k, 1) = 1/ (N + k). In 
all representations below we refer to 15 expansion coefficients a k (p), i.e. K = 14. In a few cases 
a representation with k ^ is chosen. We calculate the coefficients a k {p) in 300 logarithmic 
steps in £ e [0.4, 10 4 ] and use numerical interpolation. 

In Figures la-c we show the LO coefficient functions 4°k,' C 3°k anc ^ c fli an< ^ ^ e absolute 
accuracy of their representation by the corresponding MINIMAX-polynomial (16). Figures Id, e 
show examples for the MINIMAX-representation for NLO contributions to the heavy flavor Wilson 
coefficients. In Figures la-e below the accuracy of the representation AC(£,r)) is evaluated as 

^C^, V ) = \Cf~^,v)-m, V )\ . (18) 

Similar results as shown in Figures la-e are obtained for the other Wilson coefficients. 

Table 1 gives a complete survey on the absolute accuracy of the representations of the Wilson 
coefficients in (6,12). In many of the cases very small maximal absolute errors, i.e. of O(10~ 5 ) and 
smaller, of the Wilson coefficients were obtained in the whole £-range. In some cases the errors 
may reach ~ 6 • 10 -2 , however. Our main goal is to obtain sufficiently accurate representations 
for the heavy flavor structure functions. 
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(SU),1C = 




Figure la: LO Wilson coefficient cp L g as a function of z/p and £ = 1, 10, 10, 2 , 10 3 , 10 4 and the 
modulus of the error of the polynomial representation, AC. 
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Figure lb: LO Wilson coefficient c§l g . The MINIMAX-polynomial was determined choosing k = 0.5 
in (16). All other conditions are as in Figure la. 
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Figure lc: LO Wilson coefficient c g °^ g . All other conditions are as in Figure la. 



y 0.35 
0.3 
0.25 
0.2 
0.15 
0.1 
0.05 


^ -0.05 

S io 

| ,.; 

1 10 . 
£J 10 ' 

< IO"' 

10 " 

-I 

10 

.1 

10 

-1( 

10 



C 2)G (^TD, K = 0.4 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

z/p 



Figure Id: Contribution to the NLO Wilson coefficient Cp^ g . The MINIMAX-polynomial was determined 
choosing k = 0.4 in (16). All other conditions are as in Figure la. 
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Figure le: Contribution to the NLO Wilson coefficient dp g . All other conditions are as in Figure la. 
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Table 1: Maximal absolute errors of the MINIMAX-polynomials for the LO and NLO Wilson coefficients 
as a function of £. 
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We illustrate the principal response of the Wilson coefficients for the heavy flavor structure 
functions by the following scale-independent 'effective' distributions : 



f g (z) = 1.8z-°- 2 (l-z) 5 (19) 
f q (z) = 0.6z'°- 2 (l-z) 5 (20) 
g g (z) = 20z 1A (l-zf (21) 

j:e 2 q f q (z) = 0.15z-°- 2 (l-zf . (22) 

1,1 

for different structure functions for charm quark production choosing m c = 1.5 GeV. The strong 
coupling constant is fixed to a value of a s = 0.2 for simplicity. All convolutions are performed in 
Mellin-iV space. Finally the structure functions in x-space are obtained by the contour integral 



1 r°° 

Fi{x) = - ; 

7T JO 



/ dz\m \e^x- c ^F(c(z))] , c(z) = c + ze^ , $ ~ (3/4)tt . (23) 
Jo 1 J 



The absolute error of the MINIMAX-representation if compared to the numerical representation 
using the x-space parameterizations [14] and the sample distributions (19-22) in Figures 2a-e 
below are 

AFi(x,Q 2 ) = \F~\x,Q 2 ) - F t (x,Q 2 )\ . (24) 

As seen in the figures below the absolute errors AFi(x,Q 2 ) are mostly smaller than 10~ 5 and 
reach 10~ 4 for gi(x, Q 2 ) in a range of x < 10" 2 . 




Figure 2a: F}f% in dependence of a; and Q 2 and the absolute accuracy of the MINIMAX-representation 
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Figure 2b: F 2 L ° in dependence of x and Q 2 and the absolute accuracy of the MINIMAX-representation 
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Figure 2c: (?^o in dependence of x and Q 2 and the absolute accuracy of the MINIMAX-representation 
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Figure 2d: Contribution due to (£,r)) to in dependence of x and Q 2 and the absolute 

accuracy of the MINIMAX-representation AF£. 
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Figure 2e: Contribution due to 4 ^(£,77) to -F 2 N c LO in dependence of x and Q 2 and the absolute 
accuracy of the MINIMAX-representation AF 2 C . 



10 



4 Conclusion 



We derived semi-analytic representations for the Mellin transforms of the Wilson coefficients 
for heavy flavor production in unpolarized and polarized deeply inelastic scattering to two-loop 
order for complex values of N. These representations are obtained using Riemann-Liouville 
fractional integrals and are meromorphic functions. The expansion coefficients were determined 
using the MINIMAX-method at fixed values of £ = Q 2 /m 2 , which allows to obtain very precise 
approximations. This representation of the one- and two-loop massive Wilson coefficients is used 
in QCD evolution programs solving the renormalization group equations for mass factorization in 
Mellin— space. These evolution programs are particularly fast due to the fact that all coefficient 
functions and anomalous dimensions are only once calculated during the initialization of the code. 
In the case of the heavy flavor coefficient functions the whole array in £ is stored for the inversion 
contour as a grid allowing for precise and fast numerical interpolations during the fit of the non- 
perturbative parton distribution functions. This representation of the Wilson coefficients can 
be directly used to implement the heavy-flavor effects into the physical evolution kernels for 
factorization scheme-invariant singlet evolution. 

The parameterization can be obtained from www-zeuthen.desy.de/theory/research/ 
num . html 

Acknowledgment. This paper was supported in part by DFG Sonderforschungsbereich Tran- 
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